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Abstract. We use our non-conformal holographic bottom- up model for QCD 
described in [1 to further study the effect of the QCD trace anomaly on the energy 
loss of both light and heavy quarks in a strongly coupled plasma. We compute the 
nuclear modification factor Raa for bottom and charm quarks in an expanding plasma 
with Glauber initial conditions. We find that the maximum stopping distance of light 
quarks in a non-conformal plasma scales with the energy with a temperature (and 
energy) dependent effective power. 
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1. Introduction 

In [1] we presented a non-conformal holographic bottom-up model that captures some 
of the phenomenological properties of QCD at finite temperature. In this model, 
asymptotic freedom is replaced by conformal invariance above a certain UV scale. The 
main idea is to break the conformal invariance of A/" = 4 SYM field theory from the 
original AdS/CFT correspondence [2] in order to study effects from the QCD trace 
anomaly and its thermodynamical consequences as closely as possible. Our model is 
based on a bottom-up approach to this problem (for an alternative, top-down approach, 
see e.g. [3]) where we consider an effective, 5-dimensional gravity theory coupled to a 
scalar field The potential for the scalar field V{(j)) breaks conformal invariance in 
the infrared and it is engineered to holographically reproduce some thermodynamical 
properties of finite-temperature QCD [5J that have been computed on the lattice [6]. 
The space-time geometry is asymptotically AdS^, which is translated into conformal 
invariance of the dual field theory in the UV. 

2. Heavy quark energy loss and Raa 

A heavy quark of mass ttiq is dual to a string in the bulk which stretches from the 
bottom of a D4-brane located at a radial coordinate ~ ^/^q (in coordinates where 
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the boundary is at r = 0) to the black brane horizon at rn [Z]- Using the trailing 
string ansatz, we can study the heavy quark energy loss |8] in our model at various 
temperatures and momenta, and significant deviations from the results in the conformal 
limit have been observed around Tc 

Here we use this setup to compute the nuclear modification factor Raa of charm and 
bottom quarks (for similar calculations in a conformal plasma see [9]). We use Glauber 
initial conditions where the number density of binary collisions T^a and the participant 
nucleon density ppart are obtained from a Woods-Saxon nuclear distribution. A jet 
produced at a point x± in the transverse plane and moving in the azimuthal direction 
"sees" at time t the temperature T{x±,(l),t) oc [ppart ix± + {t — ti)e{(f))) /t]^^^, where 
ti is the initial jet production time, which we chose to be 1 fm/c. To compute Raa 
as a function of Prj, the final momentum of the quenched jet, one must know its 
corresponding initial momentum pT,i{PT,f, x±,(f)). To obtain that, one must, for a given 
x± and 0, solve the differential equation = ^ {p{t),T{x±, (j),t)), where ^ {p, T) 
is the energy loss as a (numerical) function of momentum and temperature obtained from 
our model. This equation is solved with the condition p{tf) = Ptj, where the final time 
tf is defined by a standard freeze-out condition, T{x±, (f),tf) = Tfo, where we chose 
Tfo = 150 MeV. The initial momentum is then obtained as pT,i = p{ti). Finally, one 
must average over the azimuthal directions and the transverse production points x±: 

^aaKPtj) - / d / wz^^ --r — K^) 

where Nf^in is the total number of binary collisions and -^^^ are distribution functions 
obtained from the FONLL production cross sections [lOj. 

We computed Raa for heavy quarks at RHIC and LHC energies with the initial 
temperatures at the center of the plasma of 265 MeV and 357 MeV and nucleon-nucleon 
inelastic cross-sections aNN = 4.2 and 6.3 fm^ [TT], respectively (Fig. 1, left plot). For 
bottom quarks, Raa fiattens out already at px ~ 10 GeV, which is a consequence of 
the asymptotic conformal invariance of our model. Also, note that the values of Raa 
for RHIC and LHC are very similar. In fact, even though the maximum temperature 
at LHC is higher than at RHIC, the spectral indices at RHIC are steeper than at LHC 
and, thus, the heavy quark Raa does not change much when one increases ^/s by one 
order of magnitude. 

For charm quarks Raa is relatively small, which is simply due to the fact that the 
trailing string energy loss in any holographic model decreases monotonically with mass. 
However, one should be careful when using the trailing string ansatz in this case since 
already for T ~ 290 MeV, th ~ Tm for the charm quark, which implies that the trailing 
string ansatz is not applicable in this case. 
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Figure 1. Left: Nuclear modification factor Raa at impact parameter b = 3 fm for charm and bottom 
quarks in an expanding plasma as a function of the transverse momentum px for RHIC (dashed curve) and 
LHC (solid curve). Center: Light quark stopping distance Ax as a function of energy E for different initial 
conditions (blue points) and the extracted envelope for the maximum stopping distance (red curve) . Right: 
Effective stopping power n^ff for light quarks as a function of energy E for two different temperatures. 



3. Light quark energy loss and the maximum stopping distance 

For light quarks, the bottom of the D4-brane ~ 1/ mq can approach the black brane 
event horizon th- In this case we follow [12] and study open strings with both endpoints 
on the flavor brane, which fall freely towards the black brane. These strings are dual 
to dressed qq pairs and the total spatial distance the string endpoint traverses can be 
identified with the stopping distance of light quarks in a strongly-coupled plasma. 

For a given initial energy of the string, this stopping distance will vary with the 
string initial conditions (initial spatial and velocity profile). However, it is still possible 
to obtain the maximum stopping distance for a given energy. In fact, it was shown in 
[I3| IT2] that this quantity scales as Ax^^"^. ~ E^^^ in a conformal plasma. However, in 
our model, this power gets modified and the effective power becomes temperature and 
energy-dependent, Axmax ~ E^^ff^'^'-^^ (Pig. 1, center and right plots). We see that for 
large enough energies this effective power asymptotes to the conformal value of 1/3 and 
also that the effective power is closer to the 1/3 limit the further away we are from T^, 
which are both consequences of the asymptotic conformal invariance of our model. 

Note, however, that this maximum stopping distance is not a typical stopping 
distance of light quarks; it is a rather crude quantity, which might be used as a 
phenomenological guideline, but should not be used to obtain the instantaneous energy 
loss that enters in, for example, calculations of Raa- Instantaneous energy loss, defined 
as the usual H^ worldsheet current, becomes a non-trivial quantity in this case since it 
is explicitly time-dependent and also depends on the point on the string on which one 
evaluates it. As shown in [12], at late times the instantaneous energy loss develops a 
Bragg-like peak but its magnitude depends on the point on the string where the energy 
loss is evaluated. 

Furthermore, the amount of instantaneous energy loss depends heavily on the initial 
conditions of the string (which are dual to the initial gauge field configuration associated 
with the qq pair). These two problems are the main obstacles that need to be overcome 
in order to obtain a consistent description of light quark energy loss that can be used 
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in the calculation of observables such as Raa- 

4. Conclusions, Outlook, and Acknowledgments 

In this proceedings we have extended our previous analysis of heavy quark energy loss in 
non-conformal plasmas started in [1] and computed the partonic Raa in an expanding 
plasma with Glauber initial conditions. The bottom quark Raa at RHIC decreases 
with pt and asymptotes to a value of about 0.2. Charm Raa is predicted to be heavily 
suppressed (< 0.1), but one should keep in mind that the trailing string model may not 
be applicable to describe charm quarks. According to our calculations, the heavy quark 
Raa X Pt remains practically unaltered when going from RHIC to LHC energies. 

We also showed that the typical effective power in the energy dependence of the 
maximum stopping distance of light quarks becomes energy and temperature dependent 
in a non-conformal plasma. A consistent treatment of instantaneous energy loss of light 
quarks has some practical difficulties and possible solutions are currently being explored 

We thank M. Mia, A. Buzzatti, R. Pisarski, and P. Chesler for helpful discussions. 
A.F. and M.G. acknowledge support by US-DOE Nuclear Science Grant No. DE-FG02- 
93ER40764. 
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